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In the first paper of this series, " Electrodynamics and the Gauss Linking Integral on 
the 3-sphere and in hyperbolic 3-space" we developed a steady-state version of classical 
electrodynamics in these two spaces, including explicit formulas for the vector-valued 
Green's operator, explicit formulas of Biot-Savart type for the magnetic field, and a 
corresponding Ampere's Law contained in Maxwell's equations, and then used these to 
obtain explicit integral formulas for the linking number of two disjoint closed curves. 

In this second paper, we obtain integral formulas for twisting, writhing and helicity, 
and prove the theorem link = twist + writhe on the 3-sphere and in hyperbolic 
3-space. We then use these results to derive upper bounds for the helicity of vector fields 
and lower bounds for the first eigenvalue of the curl operator on subdomains of these 
two spaces. 

An announcement of these results, and a hint of their proofs, can be found in the 



Math ArXiv, math.GT/0406276, while an expanded version of the first paper, with full 



proofs, can be found at |math.GT/ 0510388, 

The flow of this paper is indicated by the following list of sections. The first two are 
devoted to a summary of information from the preceding paper. 

1. Linking integrals in R 3 , S 3 and H 3 . 

2. Magnetic fields in R 3 , S 3 and H 3 . 

3. Link, twist and writhe in S 3 and H 3 . 

4. Proof scheme for LINK = TWIST + WRITHE. 

5. Some geometric formulas on S 3 . 

6. Proof of LINK = TWIST + WRITHE in S 3 . 

7. Proof of LINK = TWIST + WRITHE in H 3 . 

8. Helicity of vector fields on S 3 and H 3 . 

9. Upper bounds for helicity in R 3 , S 3 and H 3 . 

10. Hodge decomposition of vector fields. 

11. Spectral geometry of the curl operator in R 3 , S 3 and H 3 . 
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The integral formulas in this paper contain vectors lying in different tangent spaces; 
in non-Euclidean settings these vectors must be moved to a common location to be 
combined. 

In S* 3 regarded as the group of unit quaternions, equivalently as 577(2), the differential 
L yx -i of left translation by yx _1 moves tangent vectors from x to y. In either S 3 or 
H 3 , parallel transport P yx along the geodesic segment from x to y also does this. As a 
result, we get three versions for each of the formulas that appear in the theorems below. 

1. Linking integrals in R 3 , S 3 and H 3 . 

Let K\ = {x(s)} and K 2 = {y{t)} be disjoint oriented smooth closed curves in either 
Euclidean 3-spacc R 3 , the unit 3-sphere S* 3 , or hyperbolic 3-space H 3 , and let a(x, y) 
denote the distance from x to y. 



dy/dl 




Figure 1: Two linked curves 

Carl Friedrich Gauss, in a half-page paper dated January 22, 1833, gave an integral 
formula for the linking number in Euclidean 3-space, 

It will be convenient for us to write this as 

Lk(K 1 ,K 2 )= [ ^ x ^.V y¥ »(x,y)dsdt, 

where <p(a) = l/iAixa), and where we use </?(x, y) as an abbreviation for <^(a(x, y)). The 
subscript y in the expression V y y?(x, y) tells us that the differentiation is with respect 
to the y variable. 
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The following theorem from our first paper gives the corresponding linking integrals 
on the 3-sphere and in hyperbolic 3-space. Since the location of the tangent vectors is 
now important, we note that the vector V y y?(x, y) is located at the point y. 

Theorem 1.1. Linking integrals in S 3 and H 3 . 

(1) On S 3 in left-translation format: 

T , trs x f T dx dy If T dx. dy , , 

where (fi(a) = (n — a) cot«/(47r 2 ). 

(2) On S 3 in parallel transport format: 

Lk(K u K 2 ) = [ P y J^ x ^ • V y y?(x, y) ds dt, 
JK 1 xK 2 ds at 

where (fi(a) = (n — a) esc a/ (47T 2 ) . 

(3) On H 3 in parallel transport format: 

Ut(K u K 2 ) = [ P y J^ x % ■ V y <p(x,y)dsdt, 

JK 1 xK 2 as at 

where (p(a) = cscha/(47r). 

Greg Kuperberg (2008) obtained, independently and by a totally different argument, 
an expression equivalent to formula (2) above. 

The kernel functions used here have the following significance. 

In Gauss's linking integral, the function — tp(a) = — l/(47ra), where a is distance 
from a fixed point, is the fundamental solution of the Laplacian in R 3 , 

-Ay? = 5. 

Here S is the Dirac ^-function. 

In formula (1), the function — <f(a) = —(n — a) cot a/(47r 2 ), is the fundamental 
solution of the Laplacian on S* 3 , 



-Ay? = 5- 




Since the volume of S 3 is 27r 2 , the right-hand side has average value zero. 
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In formula (2), the function — <f(a) = — {it— a) esc a/ (in 2 ) is the fundamental solution 
of a shifted Laplacian on S 3 , 

-(A<p-<p) = 8. 

In formula (3), the function — (p(a) = — cscha/(47r) is the fundamental solution of a 
shifted Laplacian on H 3 , 

-(Ap + p) = 5. 

Our proof of the formula link = twist + writhe will depend on the asymptotic 
properties of ip at its singularity. For example, in the case of S 3 in parallel transport 
format, 

p>(a) = —^(n — a) csc(a) = — — + ——a — ——^a 2 + a 3 f (a), 

ry ' 4tt 2V ; y ' Ana An 2 2An 24tt 2 j y h 

where f(a) is bounded and smooth. Likewise 

111 7 

(p'(cx) = — - H : -a + a 2 — a 3 q(a) 

^ V ; Ana 2 2An 12vr 2 480vr yv ; 

and 

117 7 

ip (a) = ^ = H — a »a + a h(a), 

Y y ' 2na 3 Yin 2 240vr 120tt 2 v h 

where g and h are also bounded and smooth. Note that <p has no singularity at a = ir, 

in fact, (p is smooth and even around a = n: 

near ir. This implies that V y y?(a:(x, y)) exists and is zero when y is the antipodal point 
of x, even though V y a is not defined there. Because of this, the functions /(a(x, y)), 
g(a(x, y)) and /i(a(x, y)) defined above are defined, smooth and bounded for all x and 
y such that x^y. 

Because we do not need so many terms of these expansions, we will simply write: 
= + /(«), = ~ + 9(a), p"(a) = -L + h(a), 

where these new functions /, g and h are bounded and smooth everywhere on S* 3 . 
Similar calculations show that, for p(a) = cscho;/(47r) on H 3 , we again have 

= + /(«), = --^ + </(«), = ^3 + 

where these latest functions /, c/ and h are bounded and smooth everywhere on H 3 . 
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2. Magnetic fields in R 3 , S 3 and H 3 . 

In Euclidean 3-space R 3 , the classical convolution formula of Biot and Savart gives 
the magnetic field BS(v) of a compactly supported current flow v: 



BS(v)(y) = y Ei V(x)x i ^L 5d x. 

For simplicity, we write dx to mean dvol x . 

The Biot-Savart formula can also be written as 

BS(v)(y)=/ v(x) x V y y? (x,y)dx, 

7lR3 

where <po(a) = — l/(4ira) is the fundamental solution of the Laplacian in R 3 . 

In R 3 , if we start with a smooth, compactly supported current flow v, then its mag- 
netic field BS(v) is a smooth vector field (although not in general compactly supported) 
which has the following properties: 

(1) It is divergence-free, V • BS(v) = 0. 

(2) It satisfies Maxwell's equation 

V y x BS(v)(y) = v(y) + V y / v(x) • V x y? (x, y) dx, 

where <po is the fundamental solution of the Laplacian in R 3 . 

(3) BS(v)(y) ->■ as y ->■ oo. 

To see that (2) is one of Maxwell's equations, first integrate by parts to rewrite it as 

V y xBS(v)(y)=v(y)-V y / (V x • v(x)) v?o(x, y) dx. 

If we think of the vector field v(x) as a steady current, then its negative divergence, 
— V x • v(x), is the time rate of accumulation of charge at x, and hence the integral 



-V y / (V x • v(x)) v?o(x,y)dx 
Jr 3 



is the time rate of increase of the electric field E at y. Thus equation (2) is simply 
Maxwell's equation 

^ „ dE 
Vx fi = v + — . 

at 

In R 3 , S* 3 and H 3 , a linear operator satisfying conditions (1), (2) and (3) above will be 
referred to as a Biot-Savart operator. 
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Remarks. 

• To see that equation (2) above is Maxwell's equation, we integrated by parts, in 
spite of the fact that the kernel function (fo(a) has a singularity at a — 0. We 
leave it to the reader to check that the validity of this depends on the fact that the 
singularity of ip is of order 1/a. We will use this throughout the paper, without 
further mention. 

• Recall Ampere's Law: Given a divergence-free current flow, the circulation of the 
resulting magnetic field around a loop is equal to the flux of the current through 
any surface bounded by that loop. This is an immediate consequence of Maxwell's 
equation (2) above, since if the current flow v is divergence-free, this equation says 
that V x BS(v) = v. Then Ampere's Law is just the curl theorem of vector calculus. 

In particular, if the current flows along a wire loop, the circulation of the resulting 
magnetic field around a second loop disjoint from it is equal to the flux of the 
current through a cross-section of the wire loop, multiplied by the linking number 
of the two loops. Thus linking numbers are built into Ampere's Law, and once we 
have an explicit integral formula for the magnetic field due to a given current flow, 
we easily get an explicit integral formula for the linking number. 

• In R 3 , conditions (1), (2) and (3) are easily seen to characterize the Biot-Savart 
operator, as follows. Since conditions (1) and (2) specify the divergence and the curl 
of BS(v), the difference BSi(v) — BS 2 (v) between two candidates for the Biot-Savart 
operator would be divergence-free and curl-free. Since R 3 is simply connected, this 
difference would be the gradient of a harmonic function. Hence the components of 
this gradient must also be harmonic functions. Since they go to zero at infinity, 
they have to be identically zero. Thus BSi(v) = BS 2 (v). 

• In S* 3 , conditions (1) and (2) alone suffice to characterize the Biot-Savart operator, 
since there are no non-zero vector fields on S 3 which are simultaneously divergence- 
free and curl-free (i.e., there are no non-constant harmonic functions). 

• In H 3 , it is not yet clear to us how to characterize the Biot-Savart operator. Even 
strengthening (3) to require that BS(v)(y) go to zero exponentially fast at infinity 
is not quite enough. And in if 3 , unlike R 3 , the field BS(v) is not in general of class 
L 2 . 
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Theorem 2.1. Biot-Savart integrals in S 3 and H 3 . Biot-Savart operators exist 
in S 3 and H 3 , and are given by the following formulas, in which v is a smooth, compactly 
supported vector field: 

(1) On S 3 , in left-translation format: 

BS(v)(y) = / L yx -iv(x) x V y v?o(x,y)rfx- / L yx -iv(x)dx 
Js 3 47r Js 3 

+ 2V y / L yx -iv(x) • V y ^i(x,y)dx, 
Js 3 

where y?o(«) = — (tt — a) cota/(47r 2 ) and (fi(a) = —a(2n — a)/(167r 2 ). 

(2) On S 3 in parallel transport format: 

BS(v)(y) = / P y xv(x) x V y ^o(x,y)rfx, 
Js 3 

where (po(ee) = — (n — a) esc a/(47r 2 ). 

(3) On H 3 in parallel transport format: 

BS(v)(y)= / P y xv(x) x VyV?o(x,y)rfx, 

J H 3 

where (po(cx) = — cscha/(47r). 

In formula (1), the function <pi(a) = —a(2n — a)/(16n 2 ) satisfies the equation 

Ay^ =<p - [<p ], 

where [ipo] denotes the average value of tpo over S 3 . The other kernel functions already 
appeared in the linking integrals in Theorem 1.1. 

In formula (3), the magnetic field BS(v)(y) goes to zero at infinity like e~ a , where a 
is the distance from y to a fixed point in H 3 . 
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3. Link, twist and writhe in S 3 and H 3 . 

In a series of three papers (1959-1961), Georges Calugareanu defined a real- valued 
invariant of a smooth simple closed curve in R 3 by allowing the two curves in Gauss's 
linking integral to come together. In the limit, the points x(s) and y(t) now run along 
the same curve, and therefore can coincide, making Gauss's integral seem improper 
because of the |x — y| 3 in the denominator. But Calugareanu noted that in this case the 
numerator goes to zero even faster than the denominator, so that the whole integrand 
goes to zero as x and y come together, and the integral converges. In (1971), F. Brock 
Fuller called this invariant, which measures the extent to which the curve wraps and 
coils around itself, the "writhing number" : 



In those papers, Calugareanu also discovered the formula link = twist + writhe, in 
which link is the linking number of the two edges of a closed ribbon, twist measures 
the extent to which the ribbon twists around one of its edges, and writhe is the writhing 
number of that edge. 

Calugareanu proved this formula under the assumption that the simple closed curve 
K has nowhere-vanishing curvature, but the basic ideas for proving the formula without 
this assumption are already present in his papers. This can be seen in sections 6 and 
7 of this paper, where the proofs we give in S 3 and H 3 follow Calugareanu's original 
proof in R 3 , but require no curvature restriction. Nevertheless, Calugareanu's formula 
without the curvature restriction was proved by James White (1969) in his thesis, using 
a totally different approach based on ideas of William Pohl (1968a, b). 

Moving on to S 3 and H 3 , we follow Calugareanu's lead and replace the two closed 
curves K\ and K 2 in the linking integrals of Theorem 1.1 by one simple closed curve. 
Again all the integrals converge, and we use them to extend the notion of writhing 
number to these spaces. 
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Definition of the writhing integrals in S* 3 and H 3 . 

(1) On S 3 in left-translation format : 

Wt l (K) := / L yx -ip x % ■ V y <p(x,y)dsdt 
JKxK as at 

— -— [ L vx -i ^— ■ -^-ds dt, 
^ 2 Jkxk Y ds dt 

where (fi(a) = (n — a) cot a/(47r 2 ). 

(2) On S 3 in parallel transport format: 

Wtp(K) := f P y J^ x % ■ V y ^,y)dsdt, 
JKxK as at 

where (p(a) = (n — a) esc a/ (47r 2 ) . 

(3) On H 3 in parallel transport format: 

Wtp(K) := f P y J^ x ^ • VM*,y)dsdt, 

JKxK as at 

where (p(a) = cscha/(47r). 

The two versions of the writhing number on S* 3 are not the same, and one can show 
that 

Wr lW= Wr fW + 1Cngt h0fg 

The parallel transport version of writhe is more intuitively satisfying, since in this version 
the writhing number of a great circle is zero. 

We turn next to the definition of "twist" . 

Let K be a smooth simple closed curve in S 3 or H 3 , parametrized by arclength s. 
Let x(s) be a moving point along K, and T(s) = x'(s) the unit tangent vector field. 

Let v(s) be a unit normal vector field along K. Our intention is to define the (total) 
twist of v along K by a formula similar to 

Tw(v) = i-^T( S )xv( S )-v / ( S )rf S , 

the formula for twist in Euclidean 3-space. 
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T(s)xV(s) 



Figure 2: Vectors in the definition of twist 

But on S 3 there are two flavors of twist, according as v'(s) is calculated as a "left- 
invariant" derivative or as a covariant derivative. If we fall back into Euclidean mode 
and write 

, f v r v(s + h) - v(s) 
v (s) = lim - , 

h->0 h 

then the vectors v(s + h) and v(s) lie in different tangent spaces, and we must move 
them together in order to subtract. If we use left-translation in the group S 3 to move 
v(s + h) back to the tangent space at x(s) which contains v(s), then the resulting limit 
is the left-invariant derivative v' L (s). If we use parallel transport to move v(s + /i) back, 
then the resulting limit is the covariant derivative v P (s). 

The two flavors of twist on S* 3 are then given by 



Twl(v) = ± I T{s) x v(a) • v' L (s)ds (1) 



and 



One can show that 



Twp(v) = -j- I T(s) x v(s) • v'p(s) ds. (2) 
27r Jk 

length of K 



Twi(v) = Twp(v) 



2vr 



Example. Consider the great circle K = {(coss, sins, 0, 0) : < s < 2tt} on S* 3 , and 
along it the unit normal vector field v(s) = (0, 0, cos s, sin s). Then we have Tw^(v) = 
and Twp(v) = 1. 



In hyperbolic 3-space if 3 , we have only the parallel transport version of twist, 



Twp(v) = i- / T(s) x v(s) ■ v P (s) ds. (3) 
2tt ./a- 
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Now consider in S 3 or H 3 a narrow ribbon of width e obtained by starting with a 
simple closed curve K = {x(s)} and then exponentiating a unit normal vector field v(s) 
along K . One edge of this ribbon is the original curve K, and the other edge is the 
curve K e = {y £ (s)}, given explicitly (see section 5) by 

y e (s) = cose x(s) + sine v(s) in S* 3 ; 

y e (s) = coshe: x(s) + sinhe v(s) in ff 3 . 
Since if is simple, the ribbon will be embedded in S 13 or if 3 provided e is small enough. 

The vector field v(s) then points "across" the ribbon. 




Figure 3: A ribbon, its generating curve, and its vector field 
Theorem 3.1. link = twist + writhe in S* 3 and H 3 . 

(1) On S 3 in left-translation format: 

Lk(K, K e ) = Twi(v) + Wr L (K). 

(2) On S 3 in parallel transport format: 

Lk(K, K £ ) = Twp(v) + Wt p (K). 

(3) On H 3 in parallel transport format: 

Lk(K, K £ ) = Twp(v) + Wvp(K). 

We give an overview of the proof in the next section. 
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4. Proof scheme for LINK = TWIST + WRITHE. 

In spirit, our proof of Theorem 3.1 for ribbons in S* 3 and H 3 follows Calugareanu's 
original proof in R 3 : we begin with the linking integrals given in Theorem 1.1 for the 
edges K and K £ of our ribbon, let e shrink to zero, and observe the behavior of the 
linking integrand. 

The value of the linking integral is independent of e for e > since the ribbon is 
embedded and since the linking number is invariant under homotopies which keep the 
two curves disjoint. But the linking integrand blows up as one approaches the diagonal 
of K x K, and this is handled as follows. 

Outside an appropriately chosen neighborhood of the diagonal, the linking integrand 
converges to the writhing integrand as e — > 0, and its integral converges to the writhing 
number of the curve K. Inside this neighborhood of the diagonal, the linking integrand 
blows up, but its integral converges to the total twist of the normal vector field v along 



The crucial thing, recognized by Calugareanu, is that the width of the neighborhood 
of the diagonal in K x K must go to zero much more slowly than the width e of the 
ribbon. In fact, we will choose the neighborhood of the diagonal to have width e p , where 
<p < 1/3. 

To give a sense of this in action, we will outline here the proof of Theorem 3.1, part 
(2), dealing with link = twist + writhe in parallel transport format on S 3 . The 
proofs for H 3 and for left-translation format on S* 3 are essentially the same. In particular, 
in left-translation format, the integrand of the second integral in the expression for the 
linking number converges uniformly to the corresponding integrand for the writhing 
number. 

Consider, in parallel transport format on S* 3 , the linking integrand of K with K £ , 



K. 




ds 




and the writhing integrand of K, 




ds 




where <p(a) = {it — a) esc a/ (47r 2 ). 
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Then the linking number of K and K £ is given by 

Lk(K,K e )=[[ F £ (s,t)dsdt, 

J J0<s,t<L 

and the writhing number of K is given by 



Wr P (K) = [ [ F (s,t)dsdt. 

J J0<s,t<L 



Since a(x, y) is the distance between x and y, and since ip has a singularity just at 
a = 0, the only difficulty in considering the convergence of the linking integral as e — > 
happens near the diagonal, where s — t. 

To handle this, we first show that because the singularity of (p"(cx) at a = is like 
1/a 3 , we have that for sufficiently small e > 0, 



\F e (s,t) - F Q (s,t)\ <Ce 1 -^, 



provided that \s — t\ > e p . 



If < p < 1/3, then 1 - 3p > 0, and hence Ce 1 -^ as e 0. Therefore F £ {s, t) 
converges uniformly to F (s,t) in the region \s — 1\ > e p , and this region expands to the 
region \s — t\ > as e — > 0. Since the writhing integrand F (s, t) remains bounded even 
along the diagonal, this shows that 



F £ {s,t)dsdt > // F (s,t)dsdt = Wr P (K), 



|a-t|>eP J J0<s,t<L 

that is, a portion of the linking integral converges to the entire writhing integral as 
e — > 0. This is the content of Proposition 6.3 below. 

The more delicate part of the argument is the integral near the diagonal. A careful 
analysis reveals that for < p < 1, 

f t+£P 1 
lim / F £ (s, t) ds = — x'(t) x v(t) • v' P (t). 

Jt-ev 27T 

Hence 



|s-t|<eP 



1 f 

F £ {s,t)dsdt > — / x(t) x v(t) • v' P (t)dt = Tw P (v). 

2ti" Jo 



That is, the remaining portion of the linking integral converges to the entire twisting 
integral. This is the content of Proposition 6.4. In this way, we see that 

Lk(K, K £ ) = Twp(v) + Wr P (K). 
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5. Some geometric formulas on S* 3 



Before we can proceed with the details of the proof of link = twist + writhe, we 
need to collect some basic geometric formulas on S* 3 , which are treated in more detail 
in our (2008) paper. 

We consider S 3 C R 4 in the usual way, as the set 

{xG R 4 | (x,x) = 1}, 

where (x, y) is the standard inner product on R 4 . Since the linking, twisting and 
writhing integrands involve cross-products of vectors, we remind the reader that if 
xGS 3 , and v, w G T X S 3 , we define the cross product by 



V X w 



det 





Xi 


%2 


X 3 


vo 


V\ 


V2 


V3 


w 


Wi 


W 2 


W 3 


x 


Xl 


x 2 


x 3 



In this formula, we view x, v, w and the result as vectors in R 4 and {x , xi, x 2 , x 3 } 
is the canonical orthonormal basis of R 4 . From this, it is easy to see that if u is also 
tangent to S 3 at x, then the triple product u • v x w is equal to the value of the 4-by-4 
determinant whose rows are x, u, v and w. We will use the notation ||x, u, v,w|| for 
this determinant. 

Next, suppose v is a unit vector in T X S 3 . Then the unique unit-speed geodesic in S 3 
through x with initial tangent vector v is given by 

G(t) = cos t x + sin t v. 

Because (x, v) = 0, we have that (x, G(t)) = cost, and we can conclude in general that 
the geodesic distance ct(x, y) between two points x and y on S* 3 is a(x, y) = arccos(x, y). 

Moreover, if x and y are any distinct, non- antipodal points on S* 3 , then the vector 
v = (y — cos a x)/ sin a is a unit vector in T X S 3 , and the geodesic it generates connects 
x to y. From this we deduce that 

cos a x — y 



and 



V x o(x,y) 
V y o(x,y) 



sin a 



cos a y — x 
sin a 



We will also need the formula for parallel transport of a vector v G T y S 3 to T X S 3 : 

(x,v) 



V — 



l + (x,y) 



(x + y). 
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Specifically, we need the observation that P xy affects v by adding a linear combination 
of x and y. 

Using these formulas, we can make precise the definitions of F £ (s,t) and F (s,t) 
and then derive equivalent expressions for them that will be useful in our proof of 

LINK = TWIST + WRITHE. 



Proposition 5.1. Let x(s) be a simple closed curve in S 3 , let v(s) be a unit vector in 
Tx_( s )S 3 which is perpendicular to x'(s), and set y e {t) = x(i) cose + v(t) sine for each t. 
Then 



F £ {s,t) 



<ix 

ds 



xO)ye(t) 



dt 



X 



Vy,(t)¥>(a(x(s),y e (t))) 



sinov 



dy £ dx 



using the determinant notation given above and the shorthand a £ (s) = a(x(s), y £ (£)) for 
the distance between x(s) and y £ {t). Similarly, we have (for s ^ t) 



t/x / c/x 

Fo{s,t) = — -Px( s )x(t) ( x V x(t )^(a(x(s),x(t))) 



<l'(ao) 
sina 



dx 



dx 



x(t) 'dt' X ^ 



where ojq i/ie distance between x(s) and x(t). 




Figure 4: A ribbon, its generating curve, and relevant vectors 
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Proof. Using the formulas given above for V ye a and the cross product, we write 



dt 



x V ye <p(a(x,y e (t))) = (p'(a 



dy £ x cos a £ y £ - x 

dt sin a £ 

ip'(a £ ) [ dy e 

y e (t),— ,x(s) 

sin a F at 



Because the triple product is perpendicular to x and y e , this vector in R 4 is not changed 
by P X y E . Therefore we can express 



F £ (s,t) 



o?x 
ds 



sino; £ 



sin a £ 



ye(t) '^' x(s) '^ 



The proof for F is identical. 



6. Proof of LINK = TWIST + WRITHE on S 3 in parallel transport format 

In this section we prove the link = twist + writhe formula in parallel transport 
format for ribbons in the 3-sphere. As outlined above, the idea is to write the linking 
integral for the two edges of a ribbon of width e, and then take its limit as e — > 0. Of 
course the value of the linking integral stays constant, but the limit of the integral is not 
equal to the integral of the limit of the linking integrand. The latter limit is the writhe 
of the fixed edge of the ribbon, and the difference is the twist. 

To avoid unnecessary complications, we assume all our curves and deformations of 
curves are smooth, so we are free to differentiate, commute derivatives, etc. 

As indicated above, we begin with a smooth, simple closed curve K parametrized by 
arclength and given by x(s) for < s < L = length of K. We define our ribbon by 
letting v(s) be a unit vector, tangent to S* 3 at x(s) and perpendicular to T(s) = x'(s) 
for every s. The other edge of our ribbon of width e will be at distance e along the 
geodesic emanating from x in the direction of v, so it is y e (s) = cose x(s) + sine v(s). 
In general, s is not an arclength parameter for the curve y £ (s). 

The linking number of the two edges of the ribbon is: 

Lk(x,y £ )= / / F £ (s,t)dsdt, 
Jo Jo 

where F e (s,t) is given by either of the expressions in Proposition 5.1. 

The linking number is independent of e, and so our strategy will be to take the limit 
of the linking integral of the edges of the ribbon of width e as e — > 0. We will examine 
the difference between the limit of the integral (the linking number) and the integral of 
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the limit (the writhing number), and show that it is equal to the twist of the ribbon as 
defined earlier. Since the twist is defined by a single integral in contrast to the double 
integrals that define link and writhe, we'll use the following notation for the "halfway" 
integrations of the latter two quantities: 



H Lk {t;e)= [ F £ {s,t)ds 
Jo 



and 

H Wr (t) = [ F (s,t)ds. 
Jo 

Our objective will be to show that 

Mm H Lk (t; e) = H Wr (t) + "something", 

where the integral of "something" with respect to t will be the twist of the ribbon. 

As we indicated above, the convergence of the linking integrand to the writhing 
integrand fails to be uniform only near the diagonal of [0, L\ x [0, L], so we'll write 



H Lk (t) = ■■■ ds+ ■■■ ds 

J\s-t\>eP J\s-t\<eP 

where p is a number between and 1 to be determined later. We will show that the 
first term converges to H Wr (t) and the second term will give us our "something". 

Before we can prove that the convergence is uniform away from the diagonal, we need 
the following preliminary lemma. 

Lemma 6.1: There is a constant C > such that a(x(s), x(t)) > C\s — t\, where we 
consider \s — 1\ to be the "distance" on the circle with circumference L. 

Proof. This is true locally (i.e., for s near t) because x is parametrized by arclength (as 
we will justify below), and globally by compactness. 

To get the local estimate, we use Taylor's formula to write 

x(s) = x(t) + (s- t)x'(t) + ^^£ x "(0 + (s- t) 3 A(s, t) 

where A(s, t) is a bounded, smooth vector-valued function of s and t. Since x(s) lies on 
the sphere S 3 C R 4 , we have (x(s), x'(s)) = 0, and since x is parametrized by arclength, 
we have (x'(s), x'(s)) = 1. It follows that (x(s),x"(s)) = —1, and hence 

(s - 1) 2 

cosa(x(s),x(t)) = (x(s),x(t)) = 1 - v - ; + (s - t) 3 Pl (s,t), 
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where pi(s,t) is a bounded smooth scalar- valued function of s and t. In what follows, 
Pi will always stand for such a function without comment. Then clearly 

sin 2 (a) = 1 - cos 2 (a) = (s - t) 2 (l + (s - t)p 2 ) 

and using Taylor's theorem for (1 + z) l l 2 and for arcsin(,2) we conclude that 

a(x(s),x(t)) = \{s-t) + {s-t) 2 p 3 \. 

This is surely larger than ||s — t\ for \s — t\ sufficiently small, say for \s — t\ < 5. 

Corollary 6.2: a(x(s), y £ (t)) > C'\s — t\, with C independent of e, provided e is 
small enough so that the ribbon never touches itself. When \s — t\ > e p , this implies 
a(x(s),y e (t))>C'er. 

Again, this is a combination of a local estimate and a global compactness argument. 

Now we can begin to analyze the convergence of the linking integral. We start with 
the part away from the diagonal, which we expect to converge to the writhing integral. 

Proposition 6.3: IfO < p < 1/3, then 

lim II F e (s,t)dsdt= I I F (s,t) ds dt = Wr(K), 

J J \s-t\>eP JO JO 

in other words, the limit of the "away from the diagonal" part o/Lk(x, y e ) is the integral 
of Hwr{*), which is the writhing number Wr(K) . 



Proof We need to analyze the difference 

iff {a. 



F £ (s,t) - F (s,t) = 



ye(t) '^' x(s) '^ 



+ 



sinao 



x(t) '^' x(s) '^ 



sin a £ 

using the notation of Proposition 5.1. 

Using properties of the determinant, we can rewrite the difference F £ — F as a sum 



as follows: 
F £ —F = 



sina F 



dy £ dx rfx 



smceo 



dt dt 

We proceed to bound these two summands in terms of e. 

For the first summand of F s — F , 

-f'{ct £ 



^(«o) /.x tfM /,x dx 



sino; e 



dy £ dx dx 

y ^^-Tt^Ts 
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we begin with some easy preliminary observations: x(s) and dx./ds are unit vectors, and 
since y £ {t) = cose x(£) + sine v(£), we have 



dx dv 
(cose — 1 — + sine — , 
v ' dt dV 



dye _ 
dt dt 

and so can bound the second vector in the determinant as 

rfx 



dt 



dt 



< Ce, 



where C depends on the maximum value of \dv/dt\. 

To handle the first vector in the determinant, we'll group the 1/ sino; e with the y £ (t), 
and then note that the determinant is unaffected if we subtract (cos a £ / sina e )x(s) from 
the first vector. In other words, the first summand of F £ — F Q is equal to 



sin a. 



-y«(0 - 



cos a F 



-x(s), 



dy £ dx 



- ^7' x ( s )' 



dx 



sma £ dt dt ds 



And since cosa £ = (x(s), y £ (t)), the first vector in this latter determinant is a unit 
vector. Therefore, the entire determinant is bounded by \ip'(a e )\Ce. 

Finally, recall from section 1 that we can bound |v?'(oOI by a constant divided by 
a 2 , and since a £ > e p by hypothesis, we conclude that the first summand of F £ — F is 
bounded by Me l - 2p . 

The second summand of F £ — F is the determinant 

V?'(«o) (A (f/(a e ) rfx rfx 

x t - y £ (t), -,x s , — 

sin«o sinct £ dt ds 

and our job will be to handle its first vector, since the other three are all unit vectors. 

Since the value of the determinant is unaffected if we replace its first vector with 

1 , s cosa F 



<p'{a ) 



1 cosa nl , ( v 
x(t) - x(s) -<p(a e ) 



sin ao 



sin ao 



sma £ 



Ye(t)- 



sma £ 



we will obtain a bound on the determinant by bounding this vector. Using the expres- 
sions for Va derived in section 5, we write this as 

V x(s) v?(a(x(s),y e (t))) - V x(s) y?(a(x(s), x(t))) , 

and then rewrite it as 



d 

d<7 



V x(s) v?(a(x(s),y CT (t))) da, 
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where y a (t) = coscrx(t) + smav(t). 
We now calculate and estimate: 

^V x(s) v?(a(x(s),y (J (t))) = ^- (v?'(«(x(s), y CT (t))) V x(s) «(x(s), y CT (t))) 



where 
and 



da 



I = <p"(a(x(s), y ff (t))) — V x(s) a(x(s), y a (t)) 

II = <p'(a(x(s),y ff (t))) ^V xW a(x(s), y<r (t)). 



To bound |J|, we know that |Va| = 1, \da/da\ < 1, and (p"(a) — l/i^na 3 ) + 
Since we also know that a > K'e p > K'a p , we get 

Qi 



\i\ < 



a 3p ' 



To bound \II\, we have to know more about 

d , , u d / cos a CT x - y a 

— V x(s) a x (s ,y CT (t)) = t : 

da da \ sina CT 

1 da 



1 dy a cosa CT da 



sin cn^ da sin a CT rfo" sin a a da 
1 dy CT 1 da f cos a CT y CT — x 



y & 



sin a CT dcr sin dV 
1 dy a 1 dct 



sin a CT 



+ 



7 -V y<T a. 
sin a CT da sin aa 

Once again, we'll use the facts that | Va\ = 1, \da/da\ < 1, and |dy cr /(icr| = 1 to conclude 
that 



d 

da 



V x a(x,y CT ) 



< 



smct 



< 



C 



a 



Finally, we use that ip'(ct) = — l/(4na 2 ) + so that |y?'(a)l < C'/a 2 , and the fact 
that a > K'e p > K'a p to conclude that 



1 77 1 < \<p'\ 



d 

da 



Va 



< 



r 3p' 



Now we've estimated both terms into which we decomposed (d/d<7)V x ( s )</?(a:(x(s), y a (t))) , 
so we can estimate its integral as a goes from to e to obtain the result 



|V x(s) v?(a(x(s),y e (t))) - V x(s) ^(a(x(s), x(i))) | < 



Q1 + Q2 
a 3 P 



da = Qe 1 -^. 
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So far, for \s — t\ > e p , we have 

\F e (s,t) - F (s,t)\ < Me'- 2p + Qe l ~ Zp . 

Therefore 
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/ F e (s,t)ds- [ F (s,t)ds 

J\s-t\>eP J\s-t\>eP 



< L(Me x - 2v + Qe l ^ p ). 



Since 



we get 



t+sP 



F (s, i) ds 



t-eP 



< Re p , 



[ F £ (s,t)ds- [ F (s,t)ds 

J\s-t\>eP JO 

So if < p < 1/3, we can conclude that 

lim / F £ (s,t) ds = H Wr {t) 

£ ^°J\s-t\>eP 

uniformly in t, and so 



< LMe l ~ 2p + LQe 1 -^ + Re p . 



lim// F £ (s,t)dsdt = I I F (s,t) dsdt = Wt(K). 

J J \s-t\>eP JO Jo 

This completes the proof of Proposition 6.3. 

Now we must analyze the part of the linking integral near the diagonal. 
Proposition 6.4: With x 7 y £ and F £ defined as above, for < p < 1/3, 

rt+eP 



r +£ l 

lim / F £ (s, t)ds = — x'(t) x v(f) • v' P (t). 



Proof. To begin, we apply Taylor's theorem to K and write: 

x(s) = x(t) + (s - t)T(t) + ^l^!x"(t) + (s- t fA^s, t), 

where T(t) = x'(t) is the (unit) tangent vector to K at x(t) and Ai(s,t) is a smooth, 
bounded (independent of e) vector-valued function of s and t. 
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Because the link and writhe integrals (even the partial ones) are invariant under 
shifting the intervals of integration (i.e., adding different constants mod L to s and t), 
we may, without loss of generality, assume that t — 0. Then we can write: 

s 2 

x(s) = x + sT + — x" + s 3 A 1 (s) 

where x = x(0), where T = x'(0) is the unit tangent vector to K at x, where x" = x"(0), 
and where A x (s) is a smooth, bounded (independent of e and uniformly in t) vector- 
valued function of s. As in Proposition 5.1, we will write a e (s) for a(x(s), y e (0)) in what 
follows. 



Similarly, we can write 



dx. 
ds 



= T + sx" + s 2 A 2 (s), 



and we recall that the other edge of the ribbon and its derivative are given by 

y e (0) = cos e x + sin e v 

and 

cose T + sine v', 



dy £ 
dt 



t=o 



where v = v(0) and v' = v'(0). Because we are differentiating v as though it were 
a vector field in R 4 , the derivative here coincides with the covariant derivative on S* 3 , 
rather than the left-invariant derivative of section 4. Here and for the remainder of this 
section, until the statement of the theorem, we will omit the subscript in the notation 



Using the notation of Proposition 5.1, we can express 



Fe(s,0) 



sina £ (s) 

as —ip'(a £ (s))/ sina e (s) times the determinant 



cosex + sinev, coseT + sine v', x + sT + — x" + s 3 A 1 , T + sx" + s 2 A 2 
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We proceed to analyze the factor tp r (a e (s))/ sin a e (s) in front of the determinant, and 
the following four terms, into which the determinant can be expanded: 



I = 


COS £ X, 


coseT, 


6 -^" + s 3 A 1 , sx" + s 2 A 2 






II = 


COS £ X, 


sine v', 


s 2 

sT + — x" + s 3 Ai, T + sx" + s 2 A 2 




III = 


sine v, 


coseT, 


s 2 

x + — x" + s 3 Ai, sx" + s 2 A 2 






IV = 


sine v, 


sine v', 


s 2 

x + sT + — x" + s 3 Ai, T + sx" + s 2 A 2 



First, we derive an expansion of <p'(a £ (s))/ sino; e (s) in powers of s and e. To begin, 
recall that, since x(s) is a curve on S* 3 and is parametrized by arclength, we have 
(x, x) = 1, (x, T) = 0, (T, T) = 1 and (x, x") = —1 (the last equation comes from 
differentiating (x, T) = 0). Using these observations, we derive 

cosa £ (s) = (x(s),y £ (0)) 

s 2 

= (x + sT H x" + s 3 A 1 , cos e x + sin e v) 

s 2 s 2 
= cos e cos e H sin e (x", v) + s 3 p 

e 2 + s 2 

= 1 h s 3 pi + s 2 ep 2 + se 2 p 3 + e 3 p 4 

where, as before, Pi stands for a function of s and e that is bounded for all s and e, and 
smooth except perhaps for s = e = 0. 

Since sin 2 a = 1 — cos 2 a, we can conclude that 
sin 2 a £ (s) = e 2 + s 2 + s 3 p 5 + s 2 ep 6 + s£ 2 P7 + £ 3 £>8 

(S 3 S 2 £ S£ 2 e 3 \ 

1 + e^ P5 + e^T^ 6 + e^ P7 + £^ Ps ) " 

Using the Taylor series Vl + z — 1 + \z + • • • we can conclude that 

sinews) = (e 2 + s 2 ) 1 / 2 (l + ^p 9 + ^- 2 p w + -^p„ + ^Pia) ■ 
Using the Taylor series arcsin z = z + ■ • ■ we can conclude that 

/ 2 2\l/2 f-t ^ S 2 £ S£ 2 e 3 \ 
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We combine this with the expansion of (f'(a) = — 1/ (47ra 2 ) + something bounded so that 

1 



4ir(e 2 + s 2 ) 



+ P17, 



and finally conclude that 

<p'(a e (s)) ( — I 



sino; £ (s) \Att(e 2 + s 2 ) 
-1 



+ Pn 



(e 2 + s 2 y/ 2 



(1 + sp 18 + ep 19 ) 



4tt(£ 2 + s 2 ) 3 / 2 



(1 + SJ9 2 + £P2l) 



The utility of this expression for tp'(a E (s))/ sina £ (s) will become apparent when we 
multiply it by the determinants, integrate from — e p to e p , and then take the limit as 
e — > 0. Because (e 2 +s 2 ) 1 ^ 2 is larger than either s or e, we can see that whenever a+b > 3, 
the product of (p'(a £ (s))/ sino; e (s) with s a e b will integrate to something comparable to 
e p , and the integral will go to zero as e does. 

Next, we will use the observation about (e 2 + s 2 ) 1 / 2 from the preceding paragraph to 
deal with the four determinants. The first one, 



cosex, coseT, — x" + s 3 Ai, skN + s 2 A 2 



I = 

clearly has a factor of s 3 , so it will not contribute to our limit. Similarly, the second 



one, 



II 



cosex, sinev', sT H x" + s 3 A 1 , T + skN + s 2 A 2 



has a factor of s 2 sine (since you can't use the T from both the third and fourth rows), 
and so II doesn't contribute to our limit, either. 

Using the expansion cose = 1 — e 2 /2 + • • • , we can express the third determinant, 

„2 



III 

as the sum of two terms: 



sinev, coseT, x + -x" + s 3 Ai, sx" + s 2 A 5 



s sine ||v, T, x, x"|| + s 2 sine p 2 2, 



from which only the first term could contribute to our limit. 
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Finally, the fourth determinant, 



IV 



sinev, sinev', x + sT + — x" + s 3 A 1? T + sx" + s 2 A 2 



2 

can be decomposed as 

sin 2 e ||v, v',x, T|| + ssin 2 e p 2 s, 
from which only the first term could contribute to our limit. 

From our analysis so far, we conclude that 

1 



where 



Fs(s, 0) = 47r(g2 ; g2)3/2 (^ II v, T, x, x"|| + e 2 ||v, v', x, T|| + Z(s, s)), 



lim [ 5^' 9 S xl /9 ds = 0. 
J_ £V (e 2 + s 2 ) 3 / 2 



We are now ready to calculate the limit of the integral: 

lim / F £ (s,0)ds. 

From the formula for the integrand given above, this limit will equal 

( f £P se f £P e 2 \ 
lim ||v,T,x,x"|| / — — -— -ds+ ||v,v',x, T|| / — — — — - ds ). 

e^o y J„ £P Air(e 2 + s 2 ) 3 / 2 J_ £P 4vr(e 2 + s 2 ) 3 / 2 / 

The integrand in the first of these integrals is odd, so the integral is always zero (and 
hence the limit of that term is zero). For the second term, we will need the fact that 
(for p < 1) 



e 2 



Hl £P (e 2 + , 2 ) 3 / 2 dS = 2 ' 
which one calculates using the substitution x = s/e and the fact that the anti-derivative 

of 1/(1 + x 2 ) 3 / 2 is x/VTT^ 2 . 

We have thus reached our final conclusion, namely that 

lim / F e (s,t)ds = — \\vy,y:,T\\ = -(TxvV). 

£^•0 J t _ £ p lH 111 

This completes the proof of Proposition 6.4. 
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We can use Propositions 6.3 and 6.4 and a little arithmetic to start from the definitions 
of HLk(t;e) and Hy/rif) and deduce: 

Proposition 6.5: 

lim H Lk (t;e) = H Wr (t) + — (T x v v'). 

We integrate the expression in Proposition 6.5 with respect to t from to L to reach 
our final conclusion: 

Theorem 6.6: 

Lk(x,y) = Tw P (x,v)+Wr P (x). 
In other words, LINK = TWIST + WRITHE. 

Example. The simplest example of two linked curves on S 3 is a pair of great circles 
from the same Hopf fibration. We verify Theorem 6.6 in this case. The curve 
x(s) = [cos s, sin s, 0, 0] is a great circle parametrized by arclength as s runs from to 
L = 2tt. We will take x as one edge of our ribbon. 




Figure 5: The ribbon in the example 

Let v(t) = [0, 0, cost, sint]. Then v is the restriction of a left-invariant vector field to 
the great circle, and we will take the other edge of our ribbon to be 

y e {t) = cosex(t) + sine v(£) 

= [cos e cos t , cos e sin t , sin e cos t , sin e sin t] . 
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If e = 7r/2, then y £ (t) = [0, 0, cost, sin t], which is the "orthogonal" great circle to x 
and we compute the linking number of these two circles as follows. Since (x(s), y £ (t)) = 
for all s and t, we have that a(x(s), y(t)) = n/2 for all s and t. Therefore, the linking 
integrand is given by 



dx 

ds ' xye 



— — X V v ,<£ = ; 

at / sin a 



/.x rf y £ , \ ^x 

ye(t) '^' x(5) '^ 



det 



cost sint 

— sint cost 

cos s sin s 

— sin s cos s 



sin(7r/2) 
1 

The integration takes place for (s,t) G [0, 27r] x [0, 2ir], so the formula for the linking 
number of x and y e yields 1, as expected. 

To calculate the twist of our ribbon, we note that T(s) = x'(s) = [— sin s, cos s, 0, 0], 
and v'(s) = [0, 0, — sins, cos s] It is then easy to calculate that Txv-v' = [x, T, v, v'] = 1 
for all s, which gives us that the twist of the ribbon is 



Tw(x, v) 



1 

2^ Jo 



T(t) x v(t)-v'(t)dt 



1 

2^ Jo 



2n 



Idt 



To calculate the writhe of x, we use the fact that x is a geodesic, and so we have 
-fx(s)x(t)T(t) = T(s). From this it is easy to conclude that 



o?x 

ds 



x(*)x(t) 



dx. 
~dt 



X 



V xW y?(a(x( S ),x(t))) =0 



for all s and t. Therefore 



Wr(x) = 0. 



Theorem 6.6 then reads 

Lk(x, y e ) = Tw(x, v) + Wr(x) = 1 + = 1 



as it should. 
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7. Proof of LINK = TWIST + WRITHE in H 3 



The proof of link = twist + writhe in H 3 is essentially a repetition of the 
parallel transport format proof in S* 3 , except for various changes of sign and replacing 
trigonometric functions with their corresponding hyperbolic ones. In this section, we 
highlight the places where differences occur. 

As in the first paper in this series, we view H 3 C R 1,3 , the four-dimensional Minkowski 
space endowed with the inner product 

(x, y) = x y - x x y x - x 2 y 2 - x 3 y 3 

so that 

H 3 = {x G K 4 | (x, x) = 1 and x > 0}. 

We reserve the notation v • w for the induced inner product on H 3 , namely for 
v, w G T X H 3 , we define v • w = — (v,w). Because the tangent vectors are spacelike, 
this inner product provides H 3 with a Riemannian metric which is complete and has 
constant curvature —1. 

If x G H 3 , and u,v G T X H 3 , then we have 



u x v = det 



X 


X\ 


x 2 


x 3 


U 


Ui 


u 2 




v 


Vi 


v 2 




-x 


Xl 


x 2 


x 3 


U X V 


• w 


= 1 X 


u,\ 



For geodesies and the distance function, we will have 

G(t) = cosh ix + sinh t v 

for the unit-speed geodesic through x in the direction of v G T X H 3 , and the geodesic 
distance between x and y in H 3 will satisfy cosh a = (x, y ) . We have 



V y a(x,y) = 



cosh a y — x 
sinh a 



Except for the change in the inner product, the formula for parallel transport remains 
the same: the result of parallel transport in H 3 of v from y to x is 



^xy(v)=V 



(x,v) 
l + (x,y) 



(x + y). 
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Armed with these changes, and with the appropriate choice of tp(a) = csch(a:)/(47r), 
the proofs of Lemma 6.1, Corollary 6.2, Proposition 6.3 (where the biggest change is to 
have sinha e rather than sino; e in the denominator) and Proposition 6.4 proceed in the 
hyperbolic space case essentially without change from the spherical case. 

We are then led to the conclusion of Proposition 6.5, 

UmH Lk (t; e) = H Wr {t) + i(Txvv'). 

e^O ZTT 

And once again, we define the writhe of the x edge of our ribbon as 

f L f L dx f dx \ 

Wr(x) = J J ■ p x( s )x(t) ( x V x(t) ip(a(x(s),x(t))) J dsdt, 

and the twist of our ribbon as 



Tw(x, v) = — / (T x v • v') dt. 
2tt Jo 



Finally, we integrate the expressions from the hyperbolic version of Proposition 6.5 
with respect to t from to L to reach our final conclusion: 

Theorem 7.1: 

Lk(x, y e ) = Tw(x, v) + Wr(x). 
In other words, LINK = TWIST + WRITHE. 

Example. A simple example of a ribbon in H 3 C R 1 ' 3 has as one edge the circle 

x(s) = [V2, cos s, sin s, 0] 
in R 1,3 . The unit tangent vector to this curve is 

T = [0, — sin s, cos s, 0] , 
and we can choose the vector field 



y(s) = ' 



2 COS 2 S 



COS So • 

— ^,cos s, cos s sin s. sin s 



along x. Clearly, (x'(s), v(s)) = (T(s), v(s)) = for all s, and (v(s), v(s)) = — 1, so v 
is a unit vector perpendicular to T. We can make the ribbon by choosing the other edge 
to be the curve given by y e (s) = coshe x(s) + sinhe v(s). 
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Figure 6: The ribbon in the example 

By looking at the projections of the x and y £ curves into R 3 (ignoring the first 
coordinates), it's easy to see that these curves have linking number — 1. The writhing 
integrand of the x curve is easily seen to be zero, since the writhing integrand is given 
by 



x(t) 'rft' x(s) '^ 



sinha 

and the determinant is zero because the last component of each vector in the determinant 



is zero. 



For the twist of the ribbon, we must calculate T x v • v', which is given by the 
determinant 

, „ dv 



x{s),v{s), ^-,T(s) 



and calculating this determinant yields 

T x v • v' 

So we can calcuate that 



V2 



cos 2 s — 2 



Tw(x,v) 



2tt 



2tt 



(Txv- v') ds 



2- 



V2 



2n J Q cos 2 s — 2 



ds = — 1. 



Theorem 7.1 then reads 

Lk(x,y £ ) = Tw(x,v) + Wr(x) = -1 + = -1 

as it should. 
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8. Helicity of vector fields on S 3 and H 3 

Lodewijk Woltjer introduced in 1958 the notion of "helicity" of a vector field v defined 
on a domain Q in Euclidean 3-space, 



as an invariant during ideal magnetohydrodynamic evolution of plasma fields. Keith 
Moffatt (1969), recognizing that this quantity measures the extent to which the field 
lines of v wrap and coil around one another, named it "helicity" and showed that 
Woltjer's original formula could be written in the above form. 

If v is a smooth vector field on R 3 with compact support, then the above formula for 
its helicity can be written succinctly as 



where we recall that BS(v) denotes the magnetic field due to the steady current flow v. 



This is how Woltjer originally presented his invariant, J A ■ B dx, with the role of v 
played by the magnetic field B and the role of BS(v) played by its vector potential A. 

We use (8.2) to define the helicity of a vector field v on S* 3 or H 3 , and then immedi- 
ately obtain explicit integral formulas from Theorem 2.1. 

Theorem 8.3. Helicity integrals in S* 3 and H 3 . 

(1) On S 3 , in left-translation format: 




(8.1) 




(8.2) 




where <p(a) = (n — a) cot a/ (47T 2 ) and <fi(a) = — a(2ir — a)/(167r 2 ). 



(2) On S 3 in parallel transport format: 




where <p(a) = (n — a) esc a/(47r 2 ) . 
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(3) On H 3 in parallel transport format: 

H(v) = / ^yxv(x) x v(y) • V y y?(x,y) dxdy 

where (p(a) = cscha/(47r). 

In formula (1), if v is divergence-free, then the third integral in the definition of H(v) 
vanishes, and this formula then resembles the linking formula (1) of Theorem 1.1. For- 
mulas (2) and (3) already resemble the corresponding linking formulas of Theorem 1.1. 

In formulas (1) and (2), if the smooth vector field v on S* 3 is divergence-free, then 
its helicity is the same as its asymptotic (or mean) Hopf invariant, as defined by Arnold 
(1974), and is invariant under the group of volume-preserving diffeomorphisms of S* 3 . 

In formula (3), we assume that v has compact support in order to guarantee conver- 
gence of the integral. 



9. Upper bounds for helicity in R 3 , S 3 and H 3 . 

Let Q be a compact, smoothly bounded subdomain of R 3 , S 3 or H 3 , and let v be 
a smooth vector field defined on Q. Thinking of v as a current flow, its magnetic field 
BS(v) is defined by the same formulas as in Theorem 2.1, except that the integration is 
carried out only over Q. 

For uniformity of approach, we ignore the left-translation format on S 3 and write 

BS(v)(y) = / P yx v(x) x V yV ?o(x,y) dx, (9.1) 
Jn 

where 

in R 3 we have ipo(ct) = — - — so Ao? = 5; 

Ana 

in S* 3 we have (po(a) = ——^{n — a) esc a so A(p — tp = 5; 

An 1 

in H 3 we have <^o(a) = — — cscha so A(^ + y^o = 5. 

An 

The magnetic field BS(v) is defined throughout the ambient space. It is continuous 
everywhere, but its first derivatives suffer a discontinuity as one crosses the boundary 
of Q. This is a familiar situation from electrodynamics in Euclidean 3-space. 
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In what follows, we will restrict BS(v) to Q, and ignore its behavior outside this 
domain. 

Let VF(fi) denote the space of all smooth vector fields on f2, with the L 2 inner 
product 

(v, w) = / v • w d vol, 
Jn 

and associated energy (v, v) and norm |v| = (v^) 1 / 2 . 

We seek a bound for the energy or norm of the output magnetic field BS(v) on Q in 
terms of the input current flow v. Or to put it another way, we seek an upper bound 
for the L 2 -operator norm of the Biot-Savart operator, 

BS: VF(fi) ->■ VF(fi), 

in terms of the geometry of the underlying domain Q. 

As a consequence, we will determine an upper bound for the helicity H(v) = (BS(v), v) 
of the vector field v in terms of its energy (v, v) and the geometry of Q. 

Theorem 9.2. Let Q be a compact, smoothly bounded subdomain ofH 3 , S 3 or H 3 and 
let R = R(Q) be the radius of a ball in that space having the same volume as Q. Let v 
be a smooth vector field defined on Q. Then 

|BS(v)| < iV(i2)|v|, 

where 



in R 3 


we 


have 


N(R) 


= R 


in S 3 


we 


have 


N(R) 


= -(2(1 - cosR) + (ir-R) sini?) 

7T 


m H 3 


we 


have 


N(R) 


= sinh R. 



It follows immediately that the helicity H(v) = (BS(v),v) is bounded by 

|H(v)| < N(R)\v\ 2 . 

In R 3 , the overestimate N(R) = R for the norm of the Biot-Savart operator grows 
like the cube root of the volume |vri? 3 of Q. 

By contrast, in H 3 the overestimate N(R) = sinhi? for the norm of the Biot-Savart 
operator grows like the square root of the volume 27r(sinhi?coshi? — R) of Q. 
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Figure 7: N(R) for H 3 , R 3 and S 3 



Setting up for the proof of Theorem 9.2. 

To begin, let ip(a) be a real-valued function of the real variable a > 0, where we 
think of a as distance from a fixed point (and on S 3 , we have the additional condition 
< a < 7r). Assume ip has the property that 



N n (i/j) := max / |^(x,y)| d: 



is finite, where as usual we write ip(x,y) as an abbreviation for t/)(ck(x, y)). We note 
explicitly that the point y need not be chosen in fl. 

Proposition 9.3. Under the above circumstances, the operator 

T^: VF(fi) -> VF(fi) 

defined by 

T v( v )(y) = / P y X v(x) x ^(x,y)V y a(x,y)rfx 
zs a bounded operator with respect to the L 2 -norm, and 

|T^(v)|<iV n (^)|v|. 



The proof of this proposition in the R 3 case can be found in our (2001) paper, 
Lemma 3 on pages 897 and 898. The argument there follows along the lines of the usual 
Young's inequality proof that convolution operators on spaces of scalar-valued functions 
are bounded; see Folland (1995) page 9, or Zimmer (1990) Proposition B.3 on page 10. 
The proof carries over to the S 3 and if 3 cases with virtually no changes. 
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We want to apply Proposition 9.3 to the Biot-Savart operator (9.1), which we write 
BS(v)(y)= / P yx v(x)x^ (x,y)V y a(x,y)rfx = T^(v)(y), (9.4) 



where 



in R 3 we have (p' (cx) = 



Ana 2 

in S 3 we have (p' (cx) = ^(csc a + {n — a) esc a cot a) ; 

in H 3 we have (p' n (a) = — cschacotha. 

Atx 

Then by Proposition 9.3 we have 
Proposition 9.5. | BS(v)| < N n {^' )\v\. 

We turn next to estimating Nn((p f ). 
Lemma 9.6. If if? (a) is a positive, decreasing function of a, then 

Na(ip) = max / ip(yL,y)dx 
y J x en 

is maximized over all subdomains Q ofH 3 , S 3 or H 3 having a given volume when Q is 
a round ball and y is its center. 

We leave the proof of this, as well as that of the next elementary lemma, to the 
reader. 

Lemma 9.7. The functions 

V ? o(«) = 7^2 «G(0,oo) 
<p' (a) = ^— ^(csc a + (ir — a) esc a cot a) a G (0, tt] 
<p'o(a) = — csch a coth a a G (0, oo) 

47T 

are decreasing functions of a on their respective domains. 
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In view of Lemmas 9.6 and 9.7, we next compute Nn((p' ), where Q is a round ball 
of radius R in R 3 , S* 3 or H 3 , and (p' (a) is as given above. We use the shorthand 
N(R) = iW ). 

Proposition 9.8. Let VL be a round ball of radius R in R 3 , S 3 or H 3 . Then 
in R 3 we have N(R) = R 

in S 3 we have N(R) = - (2(1 - cos R) + Ok - R) sin R) 

71 

in H 3 we have N(R) = sinhi?. 

Proof. We give the proof in S* 3 and leave the other two cases to the reader. 
N(R) = / (p' (a) An sin 2 a da 

Ja=0 
f R 1 

= / — — (esc a + Ok — a) esc a cot a) An sin 2 a da 
J a =o 47r 

i r R 

— — I (sin a + Ok — a) cos a) • 

^ Ja=0 



R 

, da 

1 (2(1 - cos R) + Ok - R) sin R). 



71 



Remark. If we put R = n, then we get A^ S 3(^' ) = A/tt, in which case 

|BS(v)|<iV>(^ )|v| = -|v|, 

7T 

for smooth vector fields v defined on the entire 3-sphere. 
We contrast this with the sharp estimate 

|BS(v)| < ||v|, 

with equality if and only if v is a vector field of constant length tangent to a left or right 
Hopf fibration of S 3 . See our (2008) paper for details. 
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Proof of Theorem 9.2. By Lemma 9.7, the functions tp' (a) in R 3 , S 3 and H 3 are all 
positive and decreasing. Hence by Lemma 9.6, the quantity 



is maximized over subdomains f2 having a given volume when f2 is a round ball and y 
is its center. The values of N^(if' Q ) in that case were calculated in Proposition 9.8, and 
inserting them into the estimate 



of Proposition 9.5, we get Theorem 9.2. 

10. Hodge decomposition of vector fields 

In this section we collect, without proof, some information about the topology of 
compact subdomains in R 3 , S* 3 and H 3 , and about the structure of the space of vector 
fields on such domains. The reader will find the details in our (2002) paper. 

Let Q be a compact, smoothly bounded domain in R 3 , S 3 or H 3 , and VF(f2) the 
space of all smooth vector fields on Q, with the L? inner product and associated energy 
and norm, as defined in the preceding section. 

Let K(f2) C VF(fi) denote the subspace consisting of vector fields which are divergence- 
free and tangent to the boundary of Q, 



where n denotes the unit outward normal vector field along the boundary dfl of Q . 
These vector fields are just the incompressible fluid flows within a bounded domain, 
and in real life are naturally tangent to the boundary. In the traditional passage from 
geometric knot theory to fluid dynamics, a knot is modeled by such a flow within a 
tubular neighborhood of itself, and the flows are then called fluid knots, accounting for 
the "K" in the notation K(Q) . 

Let G(f2) C VF(f2) denote the subspace of gradient fields, 

G(f2) = {v G VF(f2) : v = Vy? for some smooth function </?: Q — > R}. 

Then we have an L 2 -orthogonal direct sum decomposition 




BS(v)|<JV n (^ )|v 



K(Q) = {v G VF(ft) : V • v = 0, v • n = 0}, 



VF(fi) = K(fi) ©G(fi). 



(10.1) 



The spaces VF(f2), K(f2) and G(f2) are all infinite-dimensional. 
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Let HK(f2) C K(f2) denote the subspace of vector fields which are not only divergence- 
free and tangent to the boundary, but also curl-free, 



We call the elements of HK(fi) harmonic knots. The subspace HK(fi) is finite-dimensional, 
and isomorphic to i^i(fi), the one-dimensional homology of f2 with real coefficients. 

The orthogonal decomposition (10.1), when further refined, yields the Hodge decom- 
position of VF(f2); see our (2002) paper for details. 

Let Q* denote the closure of the complement of Q in R 3 , S 3 or H 3 . Let g denote the 
total genus of dQ, that is, the sum of the genera of its components. Then, using real 
coefficients, H 1 (d£l) is a 2g-dimensional vector space, while Hi(Q.) and H 1 (£l*) are each 
^-dimensional, and we have the direct sum decomposition 



where the above homomorphisms are induced by the inclusions dfl C f2 and dQ C Q*. 
Let ai, a 2 , . . . , a g be a basis for ker(i?i(cftl) — > Hi(il)), and foi, & 2 , • • • , b g a basis for 



about any curve 7 in Q depends only on the homology class of 7. So we can denote this 
circulation by Circ(v, [7]). 

With this notation, the real numbers Circ(v, ai), Circ(v, a 2 ), . . . , Circ(v, a g ) are all 
zero, since the homology classes on dVL bound in Q. By contrast, 

(10.3) The real numbers Circ(v, b\), Circ(v, 6 2 ), • • • ,Circ(v,6 9 ) are in general not zero, 
and in fact define an isomorphism of HK(f2) — » H 9 . 



HK(ft) = {v G VF(ft) : V • v = 0, V x v = 0, v • n = 0}. 



(10.2) 



If v G HK(f2), then, since v is curl-free, its circulation 
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11. Spectral geometry of the curl operator in R 3 , S* 3 and H 3 

As before let f2 be a compact, smoothly bounded subdomain of R 3 , S* 3 or H 3 , and 
VF(fi) the infinite-dimensional space of smooth vector fields on f2 with the L? inner 
product. 

Now we are interested in curl eigenfields on f2, that is, vector fields vonil which 
satisfy V x v = Av for A ^ 0. In R 3 , these fields are used to model stable plasma flows; 
see our (1999) paper. 

Curl eigenfields exist for every value of A. For example, in R 3 , if 

v = sin \zi + cos Azj, 

then V x v = Av. 

We want to constrain the choice of vector fields v by interior and boundary conditions 
which guarantee that the curl operator on VF(f2) will have a discrete spectrum, while 
at the same time being reasonable for physical applications. Then we want to find a 
lower bound for the absolute values of the nonzero eigenvalues. 

To begin, we will restrict our attention to the subspace K(f2) of fluid knots, discussed 
in section 10. The vector fields in K(f2) are divergence-free and tangent to the boundary 
of Q. Since a curl eigenfield v with nonzero eigenvalue A is automatically divergence-free, 
the only real constraint here is that of tangency to the boundary. 

Let CK(f2) C K(Q) denote the subspace of vector fields whose curl lies in K(Q). Any 
eigenfield of the curl operator in K(Q) must lie in CK(f2), so restricting our attention 
to CK(f2) is no further constraint. 

Lemma 11.1. A vector field v G K(fl) lies in the subspace CK(f2) if and only if the 
circulation of v around small loops on dfl vanishes. 

Proof. The circulation of v around a small loop on dQ equals the flux of V x v through 
the small disk bounded by that loop. If this flux is zero for all such loops, then the 
normal component of V x v along <9f2 must be zero, telling us that V x v is tangent to 
<9f2, and hence that v G CK(f2). 

Remarks. 

(1) If the ciculation of v vanishes around small loops on dfl, then it also vanishes 
around homologically trivial loops there. 

(2) Any divergence-free vector field on which vanishes on <9f2 must lie in CK(f2). 
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The kernel of the map curl: CK(Q) — > K(f2) consists of vector fields on Q which are 
divergence-free, curl-free, and tangent to the boundary. These are the harmonic knots 
HK(f2) introduced in section 10. 

Since we are interested in the spectral theory of the curl operator, we would like to 
know when curl: CK(fi) — > K(f2) is self-adjoint with respect to the L 2 inner product, 
that is, when can we promise that 

(V x v,w) = (v, V x w) 

for vector fields v and w in CK(f2)? 

Lemma 11.2. Suppose that Q is simply connected, or equivalently, that all the compo- 
nents of dVt are 2-spheres. Then curl: CK(f2) — > K(Q) is self-adjoint. 

Proof. Recall the formula from vector calculus, 

V • (v x w) = (V x v) • w - v • (V x w), 
and integrate this over f2 to get 

/ ( V • (v x w) d vol = / ( V x v) • w d vol - / v • ( V x w) d vol 
Jn Jn Jn 

= (V x v, w) - (v, V x w). 

The left-hand side equals 

(v x w) • n d area, 



Ids 



I on 

and so the issue is to see when v x w has zero flux through <9f2. 

By Lemma 11.1, we know that v G CK if and only if v has zero circulation around 
small loops on dQ. Since Q is simply connected, dfl is a union of 2-spheres, and so v 
must have zero circulation around all loops on dfl. 

But this means that the restriction of v to dfl is a gradient field on that surface. So 
we write v|g^ = V/, where /: dQ — > R is some smooth function, and the gradient is 
the "surface gradient" on dfl. 

Likewise, w\qq = Vg for some smooth function g: dQ — > H. 

Now extend / and g to smooth functions F and G from Vl — > R, and consider the 
vector fields VF and VG defined on Vl. 
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Because the cross product of two gradient fields is always divergence-free, that is, 



Since V/ and Vg are, respectively, the tangential components of VF and VG along 
dfl, we can write 

VF(x) = V/(x) + a(x) n(x) and VG(x) = V#(x) + 6(x) n(x) 

for x G dQ. From this we can see that 

(VF(x) x VG(x)) • n(x) = (V/(x) x V<?(x)) • n(x) 

along dQ, and hence 

/ (v x w) • n d area = / (V/ x Vg) ■ n d area 

JdQ JdQ 

= / (VF x VG) • n,rfarea 



Thus (V x v,w) = (v, V x w) for all v and w in CK(fi), and so curl: CK(fi) ->■ K(fi) 
is self-adjoint when f2 is simply connected, completing the proof of Lemma 11.2. 

When f2 is not simply connected, the operator curl: CK(fi) — y K(Q) is not self- 
adjoint. So we seek further sensible boundary conditions which will make this operator 
self- adjoint for any domain Q. 

To this end, let v be a vector field in CK(Q). By definition of CK(f2), the curculation 
of v around all small loops on <9f2 vanishes. But then the circulation of v around any 
loop on dfl depends only on the homology class of that loop, giving us a linear map 



V • (VF x VG) = (V x VF) • VG - VF • (V x VG) = 0, 



we have 





0. 



Circ(v): H^dQ) -> R, 



from the one-dimensional real homology of dQ to the reals. 



To see where this is heading, let f2* denote the closure of the complement of in 
R 3 , S* 3 or H 3 , as in section 10, and let g be the total genus of <9f2. 
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Recall from section 10 the direct sum decomposition 

H^dQ) = ker(ifi(dft) -> ifi(ft)) + ker(// 1 (9fi) -> (10.2) 

which splits a 2g-dimensional space into two ^-dimensional summands. 

Now let AK(fi) C CK(fi) C K(Q) consist of all vector fields v in CK(fi) whose 
circulation vanishes around any loop on <9f2 which bounds in Q*. The subspace AK(fi) 
has codimension g in CK(fi). 

From (10.3), we get 

AK(Q) n HK(ft) = {0}. (11.3) 

We call AK(fi) the space of Amperian knots because, by Ampere's Law, the magnetic 
field due to a current running entirely within Q will have zero circulation around all 
loops on dfl which bound in Q* . 

We intend to show that the operator curl: AK(f2) — y K(Q) is self-adjoint, and 
proceed as follows. 

Start with a vector field v G VF(fi), and let BS(v) be the corresponding magnetic 
field defined throughout R 3 , S 3 or H 3 . Let the same symbol denote its restriction to f2, 
so that we may consider the operator BS: VF(fi) — > VF(fi). The magnetic field BS(v) 
is always divergence-free, but in general is not tangent to the boundary of Q. 

Now define BS'(v) to be the L 2 -orthogonal projection of BS(v) into K(f2). We are 
only going to apply BS' to vector fields v already in K(Q), so we regard this modified 
Biot-Savart operator as a map 

BS': K(Q) ->■ K(Q). 

We see from the orthogonal decomposition (10.1) that, for any v G K(f2), we have 
BS(v) = BS'(v) + the gradient component of BS(v). (H-4) 

Proposition 11.5. The image of the map BS': K(Q) — y K(Q) is the subspace AK(Q) 
of Amperian knots. 

Proof. Let v G K(f2), so that v is divergence-free and tangent to dQ. Then we have 
V x BS(v) = v. This follows from Maxwell's equation for subdomains of R 3 by Proposi- 
tion 1 of Cantarella, DeTurck and Gluck (2001), for subdomains of S 3 by Proposition 3.1 
of Parsley (2009), and similarly in H 3 . 
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By (11.4), we then also have V x BS'(v) = v. 

Now let 7 be a loop on dil which bounds the surface £ in Q*. Then the circulaion of 
BS'(v) around 7 equals the flux of V x BS'(v) = v through E, according to Ampere's 
Law. But the flux of v through E is zero, since v is confined to Vl. Thus BS'(v) C AK(fi), 
and we have shown that 

BS'(K(Q)) c AK(ft). 

To see the reverse inclusion, start with w e AK(f2), and let v = V x w e K(f2). We 
claim that w = BS'(v). 

To show this, first note that Vxw = v = Vx BS'(v), hence BS'(v) — w is curl-free, 
and therefore lies in HK(f2). 

Now we showed above that BS'(v) G AK(fi), and we have w e AK(Q) by hypothesis, 
so BS'(v) - w also lies in AK(fi). 

Therefore BS'(v) - w lies in AK(ft) n HK(fi) = {0}, according to (11.3), and so 
BS'(v) = w. This shows that 

AK(ft) c BS'(K(Q)), 

completing the proof. 

In the course of the proof, we actually showed a little more. 

Corollary 11.6. The maps 

BS': K(Q) -> AK(ft) and curl: AK(Q) ->■ K(Q) 
are inverses of one another. 

Proposition 11.7. The map curl: AK(fi) — > K(Q) is self- adjoint. 

Proof. Self-adjointness of this curl map is equivalent to self-adjointness of its inverse 
BS': K(f2) — > AK(fi), and this in turn is a consequence of self-adjointness of the Biot- 
Savart operator BS: VF(fi) — > VF(fi), which can be seen directly from its defining 
formula as follows. 
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Suppose that Q is a compact, smoothly bounded subdomain of R 3 , S 3 or H 3 , and 
that v, w G VF(fi). Then 

(BS(v),w>= / BS(v)(y)-w(y)dy 

= / ( / P yx v(x) x V y v?o(x,y)rfx) -w(y)rfy 
= / p yxv(x) x V y v?o(x,y) • w(y) dxdy 
= / v(x) x (-V x v?o(x,y)) • P xy w(y)cix(iy 
= / ^xyw(y) x V x <p (x, y) • v(x) dx, dy 
= / p y x w(x) x V y (p (x, y) • v(y) dx, rfy 

JSixSl 

= (v,BS(w)>, 

where we went from the third line above to the fourth by applying the parallel transport 
P xy to every term without changing the value of the integrand, from the fourth to the 
fifth by interchanging two terms and reversing the sign, from the fifth to the sixth by 
interchanging the variables x and y, and finally, comparing the sixth line to the third, 
moved on to the seventh line. This completes the proof of the proposition. 

Finally, we come to the desired result. 

Theorem 11.8. Let Q be a compact, smoothly bounded subdomain of H 3 , S 3 or H 3 
and let R = R(Q) be the radius of a ball in that space having the same volume as Q. 

Then curl: AK(f2) — > K(f2) is a self-adjoint operator, and for each v e AKQ), 

Ivl 



where 



in R 3 we have N(R) = R 

in S 3 we have N(R) = -(2(1 - cosP) + (tt - R) sinP) 



7T 



in H 3 we have N(R) = sinhP. 
In particular, if X is any curl eigenvalue on AK(Q), then |A| > 



N(R) 
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Proof. We already know from Proposition 11.7 that curl: AK(fi) — > K(fl) is self-adjoint. 

Let v G K(f2). Then, since BS'(v) is the orthogonal projection of BS(v) back into 
the subspace K(f2), we certainly have | BS'(v)| < | BS(v)|. 

But | BS(v)| < N(R)\v\ by Theorem 9.2, so that same bound applies to the modified 
Biot-Savart operator: 

|BS'(v)| < N(R)\v\. 
Now suppose that v G AK(Q). Then 



Hence 

and therefore 



V x v G K(Q) and BS'(V x v) = v. 



|v| = | BS'(V x v)| < N(R)\V x v|, 



|Vxv| ^iv^) |v| ' 



completing the proof of Theorem 11.8. 
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